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Abstract 

00 

£SJ ■ We study the spectrum of the linear operator L = -de — ede{sm6dg) subject to the periodic 

boundary conditions on 9 6 [— 7r,7r]. We prove that the operator is closed in L 2 ([—tt,tt}) with the 
domain in i?p er ([— tt, tt]) for |e| < 2, its spectrum consists of an infinite sequence of isolated eigen- 
values and the set of corresponding eigenfunctions is complete. By using numerical approximations 
of eigenvalues and eigenfunctions, we show that all eigenvalues are simple, located on the imaginary 
axis and the angle between two subsequent eigenfunctions tends to zero for larger eigenvalues. As a 
£Nj ■ result, the complete set of linearly independent eigenfunctions does not form a basis in „..([— tt, tt]). 

o ■ 
t> 
o 

^ (— I We address the Cauchy problem for the periodic heat equation 
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1 Introduction 



h = —hg — e(sin^/ie)e, t > 0, 
h(0) = h , 



XI) 



subject to the periodic boundary conditions on 9 £ [—tt,it]. This model was derived in the context of 
the dynamics of a thin viscous fluid film on the inside surface of a cylinder rotating around its axis in 
[3]. Extension of the model to the three-dimensional motion of the film was reported in [4j. 

The parameter e is small for applications in fluid dynamics [3] and our main results cover the interval 
|e| < 2 in accordance to these applications. For any e > 0, the Cauchy problem for the heat equation 
(jl.ip on the half-interval 9 € [0, tt] is generally ill-posed [13] and it is naturally to expect that the 
Cauchy problem remains ill-posed on the entire interval 9 £ [— tt,tt}. The authors of the pioneer work 
[3] used a heuristic asymptotic solution to suggest that the growth of "explosive instabilities" might 
occur in the time evolution of the Cauchy problem (jl.ip . 

Nevertheless, in a contradiction with the picture of explosive instabilities, only purely imaginary 
eigenvalues were discovered in the discrete spectrum of the associated linear operator 

acting on sufficiently smooth periodic functions f{9) on 9 £ [—tt, tt]. Various approximations of eigenval- 
ues were obtained in [3] by two asymptotic methods (expansions in powers of e and the WKB method) 
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and by three numerical methods (the Fourier series approximations, the pseudospectral method, and 
the Newton-Raphson iterations). The results of the pseudospectral method were checked indepen- 
dently in [IB] (see pp. 124-125 and 406-408). It is seen both in [3] and |16j that the level sets of the 
resolvent of (A — L)~ l form divergent curves to the left and right half-planes and, while true eigenvalues 
lie on the imaginary axis, eigenvalues of the truncated Fourier series may occur in the left and right 
half-planes of the spectral plane. This distinctive feature was interpreted in [3] towards the picture of 
growth of disturbances and the phenomenon of explosive instability. 

One more question raised in [3] was about the validity of the series of eigenfunctions associated to the 
purely imaginary eigenvalues of the operator L for e / 0. Although various initial conditions ho were 
decomposed into a finite sum of eigenfunctions and the error decreased with a larger number of terms 
in the finite sum, the authors of [3] conjectured that the convergence of the series depended on the time 
variable and "even though the series converges at t = 0, it may diverge later". This conjecture would 
imply that the eigenfunctions of L for e ^ do not form a basis of functions in the space H s ([—ir,ir]) 
with s > \ unlike the harmonics of the complex Fourier series associated with the operator L for e = 0. 

In this paper, we prove that the operator L is closed in L 2 ([— vr,7r]) with a domain in Hp er ([— 7r, 7t]) 
for |e| < 2, such that the spectrum of the eigenvalue problem 

is well-defined. Here and henceforth, we denote 

tfpUh^D = {/ € H\[-k,k\) ■ fW = /(-*)} • U-4) 

Furthermore, we prove that the residual and essential spectra of the spectral problem (|1.3p are empty 
and the eigenvalues of the discrete spectrum accumulate at infinity along the imaginary axis. We 
further prove completeness of the series of eigenfunctions associated to all eigenvalues of the discrete 
spectrum of L in Hp er ([— it, it]). Using the numerical approximations of eigenvalues and eigenfunctions 
of the spectral problem (|1.3p . we show that all eigenvalues of L are simple, located at the imaginary 
axis, and the angle between two subsequent eigenfunctions tends to zero for larger eigenvalues. As a 
result, the complete set of linearly independent eigenfunctions does not form a basis in i7p Cr ([— n, ir]) 
and hence it cannot be used directly in the studies of well-posedness of the Cauchy problem associated 
with the heat equation (11.11) . 

The paper is structured as follows. The domain, closed extensions, and the essential spectrum of 
the operator L are analyzed in Section 2. The discrete spectrum of the operator L is characterized in 
Section 3. Section 4 presents numerical approximations of eigenvalues and eigenfunctions of the spectral 
problem (jl.3p . Section 5 gives a summary and discusses the Cauchy problem for the heat equation 
(jl.ip . Appendix A reports an extension of the spectral problem (II. 3D into a self-adjoint problem in 
a weighted L 2 -space. Appendix B reformulates the eigenvalue problem (jl.3p as the resonance pole 
problem for a linear Schrodinger operator on an infinite line. These Appendices are not related to our 
main results and are left in the text for future work and references. 

2 General properties of the linear operator L 

It is obvious that the operator L is densely defined in L 2 ([— ir, ir]) on the space of smooth functions with 
periodic boundary conditions. However, the operator L is not closed in L 2 ([— it, it]) if the functions 
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are infinitely smooth. We shall therefore prove that the operator L admits a closure in L 2 ([— n,ir]) 
with a domain in n, tt]) (Lemmas Q] and [2]) . Properties of the eigenfunctions and eigenvalues of 

the spectral problem (|1.3|) in H^ cr ([— tt, tt]) are then studied (Lemmas and H]) and the absence of the 
residual and essential spectra is rigorously proved (Lemmas [5] and [6]) . 

Lemma 1 If Lf = F and F £ L 2 ([-tt,tt]) then f £ if* r ([— k, tt}) for \e\ < 2. 

Proof. For simplicity, we will set e > 0. The case e < is considered similarly. Let f(9) be in the 
domain of L on 9 £ [— tt, ir], such that Lf = F for any F £ L 2 ([—ir, ir]). By the method of variation of 
constants, we express f'{0) from the solution of the first-order ODE for f'{6): 



f{0) 



cot 1 / e (0/2) 



sin ( 



tan 1 / e (9 1 /2)F(9 1 )d9 1 + C 



00 



(2.1) 



where C is an integration constant for any 9q £ [ — tt, tt]. The function /'(#) is continuously differentiable 
at any # S (— tt, 0) U (0, 7r). Let us consider the behavior of f'{6) at the regular singular points 9q = 
and 9o = tt on 9 £ [0,tt]. A similar consideration holds at the regular singular points 9q = and 
#o = — on 9 £ [— 7r,0]. For 9q = and sufficiently small > 0, it follows from (12. lj) that 



tan 1 / £ (0 1 /2)F(^i)^ 1 



< 



tan 2 / e (^/2)^i / F 2 (9 1 )d9 1 < a\{9)9 2 ' t+1 , 



where lim ai(#) = 0. As a result, for > 0, it holds that 

where lim + a 2 (6>) / 0. If e < 2, /' G L 2 near 6» = if and only if C = 0. 
For 9q = tt and sufficiently small 7r — # > 0, it follows from (|2.ip that 



(2.2) 



J TT 



tan 1 / e {9 1 /2)F(9 l )d9 1 



< Htt - 9 1 y 2 / e d9 1 /"(vr - 9) 2 ' e tmi 2 l e {9 1 /2)F 2 {9 l )d9 l 
Je Je 

< a 2 (9)(7T - 9)- 2 / e+ \ 
where lim a\{9) = 0. As a result, for 9 < tt, it holds that 



I/' 



< 



oc x (6) 



1/e-l 



(2.3) 



where lim a 2 (#) 7^ 0. If e < 2, f £ L 2 near 9 = ir for any C / 0. Finally, we recall the Neumann- 
Poincare inequality on 9 £ [— tt, tt]: 

2 



i 2 <4vr 2 ||/'||i 2 + 



1 
2^ 



/(0)d0 



(2.4) 



Due to the estimates (|2.2h and (12. 3D . we can see that is bounded on 9 £ [— vr,7r] and so is ||/||/,2. 

Therefore, the solution f{9) of Lf = F with F £ L 2 ([-7r,7r]) lies in H 1 (\-tt if |e| < 2. □ 
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Lemma 2 The operator L admits a closure in L 2 ([— tt, tt]) for |e| < 2 mt/i Dom(L) C tt, 7r]). 

Proof. According to Lemma 1.1.2 in [6], if an operator has a non-empty spectrum in a proper subset of 
a complex plane, then it must be closed. The operator L has a non-empty spectrum in flp er ([— 7r, 7t]) 
since A = is an eigenvalue with the eigenfunction /o(#) = 1 £ i/p Cr ([— 7r, 7r]). We should show that 
there exists at least one regular point Ao £ C, such that 

V/ £ flJJI-TT, vr]) : || (L - A /)/|| L2 > *o||/IU» (2.5) 

for some fco > 0. In particular, we show that any Ao £ K is a regular point of L in .ffo C i?p Cr ([— 7r, 7r]), 
where 

tf = {/ £ F; er ([-7r,7r]) : jT /(0)d0 = oj . (2.6) 
By using straightforward computations, we obtain 

(f,Lf) = - r (1 + e cos 6) \f'\ 2 d9-e f\m9f'f"d9, (2.7) 

J — 7V J — TV 

where /) = J^ 7r g(0)f(9)d9 is a standard inner product in L 2 . If / E 7r, 7r]), then 

Re(/',/) = 0, Re(/',L/) = - ^ (l + | cos fl) |/'| 2 ^, (2.8) 
such that for any Ao € M it is true that 

|Re(/',(L-A J)/)|> (l-y) ll/'llk- 

By using the Cauchy-Schwarz inequality, we estimate the left-hand-side term from above 
|Re(/',(L-A /)/)| < |(/',(L-A /)/)| < \\f'\\ L 2\\(L - X Q I)f\\ L 2, 

such that 

\\(L-\oI)f\\ L z> (l-^)||/'|| i2 . (2.9) 

Using the Neumann-Poincare inequality (12. 4ft for any / £ Hq C -ffp Cr ([— vr, 7r]), we continue the right- 
hand-side of the inequality (|2.9p and recover the inequality (|2.5p for any Aq £ M with 



k = — ( 1 - — ) > 0. 
2vr V 2 / 

The estimate holds if |e| < 2. □ 

Remark 1 The formal adjoint of L in L 2 ([— tt,tv]) is L* = —edg (sinddg) + dg. According to Lemma 
1.2.1 in [Bj, the operator L* also admits a closure in L 2 ([— 7r,7r]) with Dom(L*) C Hp eT ([— tt, 7t]) for 
|e|<2. 

Lemma 3 Lei X be a isolated eigenvalue of the spectral problem Lf = Xf with an eigenfunction 
f e ^perQ- 71 ") 71 "])- Then > 
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(i) —A, A and —A are also eigenvalues of the spectral problem Lf = Xf with the eigenf unctions 
f(-8), f(8) andf(-e) in H^-n ,*]) . 

(it) A is also an eigenvalue of the adjoint spectral problem L* f* = Xf* with the eigenfunction f* = 
/(tt-0) inH^Q-n,*]). 

(Hi) X is a simple eigenvalue of Lf = Xf if and only if (/*, /) ^ 0. 

Proof, (i) Due to inversion 8 — ► —0, the spectral problem (|1.3f) transforms to itself with the transfor- 
mation A — > —A. Due to the complex conjugation, it transforms to itself with A — > A. (ii) Due to the 
transformation 9 —* ir — 8, the spectral problem (|1.3j) transforms to the adjoint problem L*/ = Xf 
with the same eigenvalue, (iii) The assertion follows by the Fredholm Alternative Theorem for isolated 
eigenvalues. □ 

Lemma 4 Let X be an eigenvalue of the spectral problem il.3\) with the eigenfunction f G it, tt]). 

Then 

= £ ~(77r ' !lm(A) = (77) ■ <2 ' 10) 

and Im(A) 7^ except for a simple zero eigenvalue A = 0. 

Proof. By constructing the quadratic form for / G -ffp Cr ([— tt, vr]), we obtain 

(/,£/) = e sm0\ffd8 - f ff'dO, (2.11) 

J —TT J —TT 

where the second term is purely imaginary since 



/ G F; er ([-7r )7 r]) : f f'fdO = \f(9)\ 2 \ e e =-, ~ f ff'M = -f f'fdd. (2.12) 

J —TT J —TT J —TT 

Moreover, the equality (|2.8p can be rewritten in the form 

Tm(A)(/',/) = Re(f',Lf) = ~ f_ {} + {^6) \f'(6)\ 2 d6 < - (l - ||/'||| 2 , (2.13) 

where the right-hand side is negative if |e| < 2 and f{6) is not constant on 9 G [— vr,7r]. Therefore, 
(/'>/) 7^ an d Im(A) 7^ 0. Finally, the constant eigenfunction f{9) = 1 corresponds to the eigenvalue 
A = and it is a simple eigenvalue since (/*, /) 7^ 0, where f*{9) = f(jr — 9) = 1 is an eigenfunction 
of the adjoint operator L* for the same eigenvalue A = 0. □ 

Lemma 5 The residual spectrum of the operator L in L 2 {[— vr,7r]) is empty. 

Proof. By a contradiction, assume that A belongs to the residual part of the spectrum of L such that 
Ker(L — XI) = but Range(L — XI) is not dense in L 2 {[— tt, tt]). Let g G L 2 ([—tt,tt]) be orthogonal to 
Range(L — XI), such that 

V/ G L 2 ([-vr, vr]) : = (g, (L - XI) f) = ((L* - XI)g, f). 

Therefore, (L* — XI)g = 0, that is A is an eigenvalue of L* . By Lemma[3jii), A is an eigenvalue of L and 
by Lemma [3^i), A is also an eigenvalue of L. Hence A can not be in the residual part of the spectrum 
of L. □ 
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Lemma 6 The essential spectrum of the operator L in L 2 ([— vr,7r]) is empty. 

Proof. According to Theorem 4 on p. 1438 in [8] , if L is a differential operator defined on the interval 
9 G (— vr,7r) = (— 7T, 0) U (0, 7r) and L± are restrictions of L on 6* G (— 7r,0) and G (0, 7r), then 
<7 e (L) = a e {L + ) U <7 e (L_), where cr e (L) denotes the essential spectrum of L. By the symmetry of the 
two intervals, it is sufficient to prove that the operator L + has no essential spectrum on 9 G (0, tt) 
(independently of the boundary conditions at 9 = and 9 = tt). It is also sufficient to carry out the 
proof for e > 0. 

Let us consider the spectral problem (jl.3p on 9 G [0, 7r] and use the transformation 

cos 9 = tanh t, sin 9 = sechi, t G R, 

such that G [0, 7r] is mapped to the infinite line t E R. Let f+(t) = f(9) on 9 £ [0, tt]. The function 
f+(t) satisfies the spectral problem 

- e/+(t) + /+(*) = Asecht /+(*). (2.14) 
With a transformation = e t ^ 2e g^(t), the spectral problem (|2.14p is written in the symmetric form 

- eg'l(t) + -g+(t) = Asecht g + (t). (2.15) 

Thus, our operator is extended to a symmetric operator with an exponentially decaying weight p(t) = 
sech(t). According to Corollary 3 on p. 1437 in [8], if L is a symmetric operator on an open interval 
(a, b) and Lq is a self-adjoint extension of L with respect to some boundary conditions at x = a and 
x = b, then a e (L) = a e (Lo). Here a = —oo, b = oo, and we need to show that the essential spectrum 
of the symmetric problem (|2.15p is empty in L 2 (M). This follows by Theorem 7 on p. 93 in [10]: since 
the weight function pit) of the problem —y"(t) — Xp(t)y(t) = on t G R decays faster than 1/t 2 as 
\t\ — ► oo, the spectrum of —y"(t) — Xp(t)y(t) = is purely discret^. □ 



3 The discrete spectrum of the linear operator L 

By results of Lemmas El HI [5J and El the spectral problem (|1.3p for |e| < 2 may have only two types 
of spectral data besides the simple zero eigenvalue: either pairs of purely imaginary eigenvalues or 
quartets of symmetric complex eigenvalues. We shall prove that there exists an infinite sequence of 
eigenvalues A which accumulate to infinity along the imaginary axis (Lemmas [7] and El) • We further 
prove completeness of the eigenfunctions associated to all isolated eigenvalues of the spectral problem 
(jl.3p (Theorem [1]). Finally, Theorem [2] formulates a sufficient condition for the set of eigenfunctions of 
the spectral problem (jl.3p not to form a basis in Hp Cr ([— tt, tt}). 

Lemma 7 Let < e < 2 and e/^nGl. For A G C, the spectral problem \1.3\) admits three sets of 
linearly independent solutions {fi{9), /2(G)} given by the Frobenius series 

-ti<9<tt: / 1 = l + ^ Cn r, f 2 = 9~^ (l + J2d n 9 n ) , (3.1) 

neN V neN / 

1 Although the spectral problem (|2.15p has an additional term Cy(t) with C > 0, this term only makes better the 
inequality (30) on p. 93 in the proof of Theorem 7 of [10] , 



6 



0<±9<n: /± = l + ^a±(7TT^r, ft = (* T U + £ h t ^ T 0) n , (3.2) 

neN V neN / 

where all coefficients are uniquely defined. The solution fi(8) is an analytic function of X G C uniformly 
on 6 € [— TT, 7r]. 

Proof. Existence of two linearly independent solutions on — 7r < 9 < tt in the form (|3.ip and on 
< ±0 < 7r in the form (|3,2p follows by the ODE analysis near the regular singular points [5j. The 
difference between the two indices of the indicial equation is | and it is non- integer for 6/i,n6 pH. 
Since the spectral problem (jl.3p depends analytically on A and the Frobenius series converges absolutely 
and uniformly in between two regular singular points, the solution fi{9) is analytic in A £ C for any 
fixed 6 G (— 7r,7r). Due to uniqueness of the solutions of the ODE (|1.3p . the solution f\{9) can be 
equivalently represented by the other solutions 

h(0) = A ± ff{6) + B ± f^(9), < ±9 < 7T, (3.3) 

where A 1 * 1 and -B 1 * 1 are some constants, while the functions ft {9) and ft (9) are analytic in A G C for 
any fixed ±# € (0, tt]. By matching analytic solutions for any ±9 G (0, tt), we find that A 1 * 1 and I? 1 * 1 
are analytic functions of A G C, the Frobenius series for /i(#) converges absolutely and uniformly on 
G [— 7r,7r], and the solution is an analytic function in A G C uniformly on 9 G [— w, tt]. □ 

Corollary 1 There exists an analytic function F e (X) on ImA > 0, roots of which give isolated eigen- 
values of the spectral problem il.3\) with the account of their multiplicity. The only accumulation point 
of isolated eigenvalues in the X-plane may occur at infinity. 

Proof. The function / G H 1 ([—tt,it]) satisfies the spectral problem (|1.3|) if and only if f(9) = Cofi(0) 
on 9 £ [— tt, it], where Cq = 1 thanks to the scaling invariance of homogeneous equations. By using 
the representation (13. 3p . we can find that = lim f\{9) are uniquely defined analytic functions in 

A G C. The function F e (X) = A + — A~ is analytic function of A G C by construction and zeros of F e (X) 
on ImA > coincide with the eigenvalues A of the spectral problem (II. 3p with the account of their 
multiplicity. If F e (Ao) = for some Ao G C, the corresponding eigenfunction f{9) lies in Hp ET ([— tt, it]), 
i.e. it satisfies the periodic boundary conditions f(ir) = /(— tt). By analytic function theory, the 
sequence of roots of F e (X) can not accumulate at a finite point on A G C. □ 

Remark 2 We will use the method involving the analytic function -F e (A) on A G C for a numerical 
shooting method which enables us to approximate eigenvalues of the spectral problem (jl.3p . This 
method involves less computations than the shooting method described in Appendix C of [3]. Nev- 
ertheless, it is essentially the same shooting method and it uses the ODE analysis near the regular 
singular point (Lemma [7]) , which repeats the arguments in Appendix B of [3] . 

Lemma 8 Fix < e < 2 and let {A n } n6 N be a set of eigenvalues of the spectral problem U.3\) with 
ImA n > 0, ordered in the ascending order of \X n \. There exists a finite number N > 1, such that for all 
n > N, X n = iuj n G iK+ and 

uj n = Cn 2 + o(n 2 ) as n — > oo, (3.4) 

for some C > 0. 

2 An additional logarithmic term log(-7r — 9) may need to be included into the Frobeneus series if e = ^, n £ Z. 
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Proof. We reduce the spectral problem (|1.3|) to two uncoupled Schroddinger equations on an infinite 
line. Let f{6) be represented on two intervals ±9 £ [0, ir] by using the transformations 

cos# = tanhi, s'm9 = isechi, (3-5) 

where t £ R. Then, the functions f±(t) = f(6) on ±8 £ [0, 7r] satisfy the uncoupled spectral problems 

- e/£(i) + = iAsecbi /±(i), t £ K, (3.6) 

The normalization condition /(0) = 1 is equivalent to the condition lim /±(i) = 1. The periodic 

t— »oo 

boundary condition /(vr) = /(— vr) is equivalent to the condition lim f—(t) = lim f+(t). The linear 
problems (|3.6p are reformulated as the quadratic Ricatti equations by using the new variables 

f ± (t) = e& S ±(*'K ■ S± - e{S' ± + S 2 ± ) = iAsechi. (3.7) 

We choose a negative root of the quadratic equation in the form 

5 ± (t)= 1 -V /lT4eAsecht - 4£2 ^, R± = S> ± {t). (3.8) 

The representation (|3.8|) becomes the chain fraction if the derivative of S±(t) is defined recursively 
from the same expression f|3 . 8|) . By using the theory of chain fractions, we claim that R± = 0(-v/[X|) 
as |A| — * oo uniformly on t £ R. The function F £ (X) of Corollary [T] is now expressed by 



F € (X) = lim [f + (t) - /_(*)] = eJ-«. 6 +W* - e J-«,f-W««. (3.9) 

t— >— oo 



J? OB s + (t)dt_j? OB s-(t)dt 

Zeros of -F e (A) are equivalent to zeros of the infinite set of functions 
1 f°° r 

GJX) = - / Jl + 4eAsecht - 4e 2 J R_(t) - Jl - 4eAsecht - 4e 2 R + (t) alt - n, (3.10) 

47rie J_ 00 L J 

where n £ N. If R±(t) = 0, the function G n (uj) = G(iuS), n £ N is real-valued and strictly increasing 
on a; £ R + with G n (0) = —n. By performing asymptotic analysis, we compute that 



I poo r 

— / J I + 4iewsecht - 4e 2 R-(t) - J I - 4ieusecht - 4e 2 R+(t) 



47rze 

1 Z" 00 2itosecht + e(i?+ - i?_ 



dt 



ni J-oo y/l + 4iewsech£ - 4e 2 _R_ (t) + ^1 - 4ieu;sechi - 4e 2 i?+(t) 
/u f°° dt 



dt 



\f2eir J-oo V cosh t ' ^ ^ 

such that lim G n (u) = oo. Therefore, there exists exactly one root to = uj n of G n {u)) for each n. Since 

i?_ = -R+ for A = iuj £ iR, each simple root of G n (uj) persists for non-zero values of R±{t) = 0(y / cJ) 
uniformly on t £ R as to — > oo. According to the asymptotic result (|3.1ip . the roots uj n of G n (u) satisfy 
the asymptotic distribution (13.41) with C = -. — r-j. □ 

V" 00 Vcosh i / 

Remark 3 Analysis of Lemma [5] extends the formal WKB approach proposed in Section 3 of [3] . In 
particular, the equation (|3.10p with R± = has been obtained in Eq. (3.11) of [3]. 
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Theorem 1 Let {f n (9)}n<=N be the set of eigenf unctions corresponding to the set of eigenvalues {X n } n &N 
in Lemma\E with ImA n > 0. The set of eigenf unctions is complete in a subspace of Hq C ilp er ([— tt, 7r]). 



Proof. By Corollary [IJ eigenvalues of L with ImA > accumulate to infinity, such that the matrix 
operator M = L _1 defined on Ho C -£fp Cr ([— tt, tt]) is compact. By LemmaEl there are infinitely many 
isolated eigenvalues and large eigenvalues are all purely imaginary, such that |A n | = 0(n 2 ) as n — > oo. 
These two facts satisfy two sufficient conditions of the Lidskii's Completeness Theorem. According 
to Theorem 6.1 on p. 302 in [11], the set of eigenvectors and generalized eigenvectors of a compact 
operator M in a Hilbert space H is complete if there exists p > such that 



-i 



s n (M) = o(n p ), as n — > oo, (3-12) 
where s n is a singular number of the operator M, and the set 

W M = {(Mf,f): feH, 11/11^2 = 1} (3.13) 
lies in a closed angle 9 m with vertex at and opening ^. 

Since the singular numbers s n are eigenvalues of the positive self-adjoint operator (MM*) 1 / 2 and the 
eigenvalues of L grow like 0(n 2 ) as n — > oo, we have s n (M) = 0(n -2 ) as n — > oo, such that the first 
condition (|3. 12|) is verified with p = 1. Since all ImA rt > for the set of eigenvalues {A n } n£ N of Lemma 
El the spectrum of M lies in the lower half plane, such that the second condition (|3.13|) is also verified 
with p = 1 (6m = I") . □ 

Corollary 2 The set of eigenf unctions {f n (9)} n ez with fo = l and /_ n = f n , Vre G N is complete in 

Remark 4 Due to linear independence of eigenfunctions for distinct eigenvalues, the set of eigen- 
functions {/n(#)}nez is also minimal if all eigenvalues are simpl^l. If the set {/ n (#)}nez is com- 
plete and minimal, any function / G H^ CT ([— tt, tt]) can be approximated by a finite linear combi- 

nation /jy (6*) = 2^ c nfn(9), such that for any e > 0, there exists N > 1 and the set of coef- 

n=-N 

ficients {c n }-N<n<N, such that the inequality ||/ — fN\\ii^ et ([— k,%]) < 6 holds. This approximation 
does not imply that the set {f n {9)} n & forms a Schauder basis in #*([— it, tt]), in which case there 
would exist a unique series representation f(6) = Yl c n f n (9) for any / £ ffp er ([— 7r, 7r]), such that 

iV lim o ll/-/iv||Hi ef ([-^ 7 r])=0. 

Theorem 2 Let {/n(#)}nez ^ e a complete and minimal set of eigenfunctions of the spectral problem 
\1.°J\) for the set of eigenvalues {X n }nez in Theorem^ The set of eigenfunctions does not form a basis 
in Hilbert space H^d— tt, tt]) if lim cos(/ n ,/ n+ i) = 1. 

Proof. The Banach Theorem defines a condition that the complete and minimal set of eigenfunctions 
{fn{9)} n £z forms a basis in Hilbert space Hp CT ([— tt, tt]). According to the Theorem 2 on page 31 in 



3 By Lemma [8] all eigenvalues are simple starting with some n > N. 
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the complete and minimal set of eigenfunctions forms a basis if and only if sup ||Pjv|| < 00 > where 

N 

P/v is the projector of the linear span {f n }-N<n<N in the direction of the linear span {f n }\ n \>N+i- 

Since the Hilbert space Hp CT ({— n, ir]) is a direct sum of the two linear spans above, the norm of the 
parallel projector Pjy has the geometrical representation ||P/v|| = 1 where a at is the angle between 
the two linear spans p]. This implies that the set {fn(6)}nez is a basis in the Hilbert space H^ cr ([— ir, ir]) 
if 

ft ~r " \ K/nj/n+l) n / Q -.a 

cos(/n,/n+i) = n , mi , i| < 1, (3.14) 

||/n||||/n+l|| 

for sufficiently large n G Z [12J. In the other words, the angles between the eigenfunctions should 
be uniformly bounded away from zero as n — > 00. If the angle tends to zero as n — > 00, the set of 
eigenvectors is not a basis in the Hilbert space. □ 



4 Numerical approximations 

We approximate isolated eigenvalues of the spectral problem (j 1 . 3 1) for < e < 2 numerically. In 
agreement with numerical results in [3], we show that all eigenvalues in the set {\ n }nez are simple and 
purely imaginary. Therefore, the set {A n } ng z can be ordered in the ascending order, such that Ao = 0, 
A n = — A_ n , Vn E N, ImA n < ImA n+ i and lim |A n | = 00. We also show that the angle between two 

subsequent eigenfunctions f n {9) and f n+ \{9) in the set {/n(#)}nez tends to zero as n — ► 00. 



4.1 Shooting method 

The numerical shooting method is based on the ODE formulation of the spectral problem fjl .3[> . By 
Lemma [7] and Corollary [lj complex eigenvalues A £ C are determined by roots of the analytic function 
F e (\) in the A-plane. The number of complex eigenvalues can be computed with the winding number 
theory. The number and location of purely imaginary eigenvalues can be found from real-valued roots 
of a scalar real-valued function. 



Proposition 1 Let the eigenfunction f(9) of the spectral problem \1.3\) for < e < 2 be normalized 
by the condition /(0) = 1. The eigenvalue A is purely imaginary if and only if f(9) = f(—9) on 
9 G [-vr,7r]. 



Proof. If A € iM and f(9) satisfies the second-order ODE (|1.3j) on 9 £ [— tt, 7r], then f{—9) satisfies the 
same ODE ((LHJ) on 9 G [-vr,vr]. By Corollary QJ if / 6 ^ er ([-vr, vr]), /(0) = 1 and < e < 2, the 
solution f(9) is uniquely defined. By uniqueness of solutions, f(9) = f(—9) on 9 £ [—tt,tt]. 

If f(9) = J{-9) on 9 e [-tt, tt], then, 

sin0|/'(0)| 2 d0= / sin0|/'(0)| 2 ^- / sin 9\ f ' {-9)\ 2 d9 = 0, 
such that ReA = according to the equality (|2.10p in Lemma 0] □ 
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Corollary 3 Let f{6) be an eigenfunction of the spectral problem 111.3]) for A £ iM, such that f £ 
ilp er ([— 7r, 7r]) and /(0) = 1. 27ien ; /(7r) = /(— 7r) is equivalent to f(n) £ R. T/ie eigenvalue X £ iR is 
simple if and only if 

(/* , /) = 2Re f /(0)/(tt - 0)d0 / 0. (4.1) 

JO 



Proof. The first assertion follows by the symmetry relation /(#) = f(—9) evaluated at 9 = ir. The 

second asserion follows by Lemma[3]with the use of the symmetry f*(9) = f(ir — 9). □ 

Numerical Method: By using Lemma [71 the function f{9) with /(0) = 1 is represented uniquely by 
the Probenius series 

f(9) = f 1 (9) = l + J2cn9 n , (4.2) 

neN 

where the coefficients {c n } ng N are uniquely defined by the recursion relation 



(-l)~m \ 

Ac n _i + en > T~nT Cm ' n G N ' ( 4 - 3 ) 

(n — m + 1 ! / 

where Co = 1 and N' is a set of integers in the interval [1, n — 2] such that n — m is even. For instance, 

A A 2 A(A 2 + e(l + 2e)) 



1 + e' z 2(l + e)(l + 2e)' ° 3!(1 + e)(l + 2e)(l + 3e) ' 

and so on. We truncate the power series expansion on N = 100 terms and approximate the initial 
value [f(9o),f'(9o)] at #o = 10~ 8 . By using the fourth-order Runge-Kutta ODE solver with time step 
h = 10 -4 , we obtain a numerical approximation of / = f+(9) on 9 £ [9o,ir — 9q] for A and / = f-(9) 
on the same interval for —A. By Lemma [3^i), the numerical approximation of the function F t (X) of 
Corollary [T] is 

F e (A) = / + (7T-0 o )-/-(7r-0o). (4.4) 

If A £ iR, the function F £ (X) is simplified by using Corollary [3] as F e (X) = 2ilmf + (-K — 9q). Table 1 rep- 
resents the numerical approximations of the first four non-zero eigenvalues A £ iM. for e = 0.5, 1.0, 1.^3 
with the error computed from the residual 



R 



(f,Lf) x 



(/,/) 



We can see from Table 1 that the accuracy drops with larger values of e and for larger eigenvalues, but 
the eigenvalues persist inside the interval |e| < 2. 

Figure 1 shows the profiles of eigenfunctions f{9) on 9 £ [0, tt] for the first two eigenvalues A = iuJip £ 
zM_l_ for e = 0.5 (left) and e = 1.5 (right). We can see from Fig. 1 that the derivative of f{9) becomes 
singular as 9 — > ir~ for e > 1. We can also see that the real part of the eigenfunction f(ff) has one 
zero on 9 £ (0, it) for the first eigenvalue and two zeros for the second eigenvalue, while the imaginary 
part of the eigenfunction f{9) has a fewer number of zeros by one. The numerical approximations of 



4 We note that the Frobenius series (|4.2|l is not affected by the logarithmic terms for e = 0.5 and e = 1.0, since is the 
largest index of the indicial equation at 8 = 0. 
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Figure 1: The real part (blue) and imaginary part (green) of the eigenfunction f(9) on 9 € [0,tt] for 
the first (solid) and second (dashed) eigenvalues A = iui^ € 2R+ for e = 0.5 (left) and e = 1.5 (right). 

the eigenvalue and eigenfunctions of the spectral problem (|1.3p are structurally stable with respect to 
variations in 9q, N and h. 

Figure 2 shows the complex plane of w = F e (\) (left) and the argument of w (right) when A traverses 
along the first quadrant of the complex plane A G Ai U A2 U A3 for e = 0.5. Here Ai = x + ir with 
x € [r, R], A2 = Re ltp with tp £ [tpo, | — ipo] and A3 = r + iy with y £ [r, R], where r = 0.1, R = 10, and 
ipo = arctan(r/i?). It is obvious that the winding number of F e (X) across the closed contour is zero. 
Therefore, no zeros of F e (X) occurs in the first quadrant of the complex plane A G C. The numerical 
result is structurally stable with respect to variations in r, R and e. 



e 


UJl 


Ri 


u>2 


R2 


0J 3 


R3 




i?4 


0.5 


1.167342 


0.000051 


2.968852 


0.000405 


5.483680 


0.001436 


8.715534 


0.003653 


1.0 


1.449323 


0.000837 


4.319645 


0.007069 


8.631474 


0.024964 


14.382886 


0.061881 


1.5 


1.757278 


0.002691 


5.719671 


0.018412 


11.846709 


0.054271 


20.138824 


0.113834 



Table 1: Numerical approximations of the first four eigenvalues A = iiv n of the spectral problem 
(jl.3p and the residuals R = R n for three values of e. 



4.2 Spectral method 

The numerical spectral method is based on the reformulation of the second-order ODE (jl.3p as the 
second-order difference equation and on the subsequent truncation of the difference eigenvalue problem. 
It is found in [16] that the truncation procedure lead to spurious complex eigenvalues which bifurcate 
off the imaginary axis. 

Numerical method: Let f{9) be an eigenfunction of the spectral problem (]1.3p in H^ CT ([— n, it]). 
This eigenfunction is equivalently represented by the Fourier series 

M = £ fne~ ine , f n = ±- r f(9)e in0 d9, (4.5) 
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Figure 2: The image of the curve w = F e (X), when A traverses along the contours Ai (blue), A2 (green) 
and A3 (magenta) for e = 0.5: the image curve on the tu-plane (left) and the argument of w (right). 



where the infinite-dimensional vector f = (..., /_2, f-i, fo, /1, f'2, •••) is defined in f £ equipped 
with the norm ||f || 2 2 = SnezU + n2 )\fn\ 2 < 00 • The spectral problem (jl.3p for |e| < 2 is equivalent to 
the difference eigenvalue problem 



nf n + -n [{n + l)/ n +i - (ra - l)/ n -i)] = -i\f n 



n S 



The difference eigenvalue problem (|4.6p splits into three parts 

Af+ = -i\f+, Af_=fAf_, A/ = 0, 
where i± = (/±i, f±2, •••) and ^4 is an infinite-dimensional matrix 



A 



(4.6) 



(4.7) 



1 


e 








— e 


2 


3e 








-3e 


3 


6e 








-6e 


4 



(4- 



Since A = D — iS, where D is a diagonal matrix and 5 is a self-adjoint tri-diagonal matrix, one can 
define the discrete counterpart of Lemma H] 



ImA 



(f+,£f+) _ Ener^l/nl 2 



ReA 



(f+,gf4 
(f+,f+ 



where ImA > 0. The adjoint eigenfunction f*{6) = /(tt — 9) is recovered from the eigenvector f by 
f* = Jf , where 

"0 J " 
J= 1 
_J 

and Jq is a diagonal operator with entries (—1, 1, —1, 1, ...). 
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Figure 3: The distance between eigenvalues computed by the shooting and spectral methods for e = 0.1. 

According to Theoremdl rewritten from the set of eigenfunctions {/ n (^)} n ez to the set of eigenvectors 
{f n }n&, the inverse matrix operator A -1 is of the Hilbert-Schmidt type. Let A N denote the truncation 
of the matrix operator A~ 1 at the first N rows and columns. If a sequence of truncated operators Aj^ 
converges uniformly to the limiting compact operator A^ 1 as N — > oo, then the spectra of matrices 
Ajf 1 also converge to the spectrum of A' 1 as N — > oo. However, the distance between eigenvalues of 
A N l and A' 1 may not be small for fixed N. 

The smallest eigenvalues of the truncated matrix A N are found with the parallel Krylov subspace 
iteration algorithm [9] . Figure [3] shows the distance between eigenvalues of the shooting method and 
eigenvalues of the Krylov spectral method for e = 0.1. The difference between two eigenvalues is small 
of the order O(10~ 3 ) but the advantage of the parallel algorithm is that the calculating time of 20 
largest eigenvalues of A N l for N = 10 6 takes less than one minute on a network of 16 processors while 
finding the same set of eigenvalues by the shooting method with the time step h = 10 -5 takes about 
one hour. 

Figure U] shows symmetric pairs of eigenvalues of the matrix An for e = 0.3 at N = 128 (left) and 
N = 1024 (right). We confirm the numerical result of [16] that the truncation of the matrix operator 
A always produces splitting of large eigenvalues off the imaginary axis. Moreover, starting with some 
number n, the eigenvalues of An are real- valued. This feature is an artifact of the truncation, which 
contradicts to Lemmas [J] and [8] as well as to results of the shooting method. However, the larger is N, 
the more eigenvalues remain on the purely imaginary axis. Therefore, the corresponding eigenvectors 
can be used to compute the angle in Theorem [2j 

Figure [5] (left) show the values of the cosine of the angle (|3.14p for the first 20 purely imaginary 
eigenvalues for e = 0.1. As we can see from the figure, the angle between two eigenvectors tends to 
zero for larger eigenvalues up to the numerical accuracy. Figure [5] (right) and Table 2 show that the 
angle drops to zero faster with larger values of the parameter e. 




14 



Numerical eigenvalues of L_ e = 0.3; N = 1024 Numerical eigenvalues of L e = 0.3; N = 128 

60 1 1 1 1 1 1 1 1 15, 1 1 1 ' 1 1 




500 1000 1500 



Figure 4: Spectrum of the truncated difference eigenvalue problem (|4.6p for e = 0.3: N = 128 (left) 
and N = 1024 (right). 




e 



Figure 5: Left: the values of cos(/„, f n +i) for the first 20 purely imaginary eigenvalues for e = 0.1. 
Right: the values of cos(/i,/2) versus e. 
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eigenvectors 



e = 0.1 



e = 0.3 



e = 0.5 



1- 2 

2- 3 

3- 4 

4- 5 

5- 6 

6- 7 

7- 8 

8- 9 

9- 10 

10- 11 

11- 12 

12- 13 

13- 14 

14- 15 

15- 16 

16- 17 






















.120166 
.461330 
.680709 
.799235 
.858944 
.892869 
.914745 
.930023 
.941262 
.949843 
.956580 
.961987 
.966407 
.970073 
.973153 
.975764 






















.325116 
.716192 
.838889 
.890440 
.921498 
.940395 
.953124 
.962120 
.968732 
.973741 
.977629 
.980702 
.983297 
.983459 
.995335 
.998749 






















.431987 
.780641 
.878055 
.914622 
.940306 
.955239 
.965235 
.972204 
.977265 
.981057 
.983988 
.986072 
.989617 
.990547 
.999101 
.999601 



Table 2: Numerical values of cos(/ n , f n +i) f° r the first 16 purely imaginary eigenvalues for three 
values of e. 

The angle between two subsequent eigenvectors is closely related to the condition number 



which measures the norm of the spectral projections [15] . By Lemma [3]^ hi) , the condition number is 
infinite for multiple eigenvalues since (f n , /*) = 0. From the point of numerical accuracy, the larger is 
the condition number, the poorer is the structural stability of the numerically obtained eigenvalues to 
the truncation and round-off errors. 

Figure [6] shows the condition number (j4.9|) computed for the first 40 purely imaginary eigenvalues 
for e = 0.001 and e = 0.002. We can see that the condition number grows for larger eigenvalues which 
indicate their structural instability. Indeed, starting with some number n, all eigenvalues are no longer 
purely imaginary. The condition numbers become extremely large with larger values of e. 

We finally illustrate that all true eigenvalues of the spectral problem ([1.3]) are purely imaginary and 
simple. To do so, we construct numerically the sign-definite imaginary type function and obtain the 
interlacing property of eigenvalues of the spectral problem (II. 3D for two values e = eo and e = e%, 
where |ei — eo| is small. We say that the eigenvalues exhibit the interlacing property if there exists an 
eigenvalue for e = e\ between each pair of eigenvalues for e = eo and vice verse. 

A meromorphic function G{X) is called a sign-definite imaginary type function if ImG(A) < 

(ImGr(A) > 0) on Im(A) < (Im(A) > 0) [2j. We construct the meromorphic function G(u) in the form 
F (X) 

G(X) = p° , where F e (X) is an analytical function of Corollary [TJ The numerical approximation of 
the meromorphic function G(X) is given by G(\) = ~ . According to Theorems II. 2.1 - II. 3.1 on p. 
437-439 in [2], the function G(A) is a meromorphic function of sign-definite imaginary type if and only 




(4.9) 
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Figure 6: The condition number for the first 40 purely imaginary eigenvalues for e = 0.001 (red) and 
e = 0.002 (blue). 

if it has the form G(X) = ^|^y where P(\) and Q(X) are polynomials with real coefficients, with real 
and simple zeros, which are interlacing. 

Table 3 shows this interlacing property of eigenvalues for eo = 0.48 and e± = 0.5. The remainder 

\\L f—\f\\ 

term R e = 1 |p^|| measures the numerical error of computations. We have also computed numerically 

the values of G(X) on the grid 0.1 < ImA < 100 and 0.1 < ReA < 100 with step size 0.1 in both 
directions (not shown). Based on the numerical data, we have confirmed that the function G(X) does 
indeed belongs to the class of sign-definite imaginary type functions while the eigenvalues {A n } ne z 
exhibit the interlacing property. This computation gives a numerical verification that all eigenvalues 
of the spectral problem (|1.3j) are simple and purely imaginary. 



ImA eo 




ImA ei 


R ei 


1.063112 


2.32438e 


- 10 


1.068314 


2.40727e 


- 10 


2.970880 


2.19667e 


- 10 


3.024428 


2.25307e 


- 10 


5.414789 


2.20243e 


- 10 


5.542829 


2.26833e 


- 10 


8.471510 


2.0904e • 


- 10 


8.693066 


2.15721e 


- 10 


12.312548 


2.00793e 


- 10 


12.665485 


2.06007e 


- 10 


16.816692 


1.97653e 


- 10 


17.327038 


2.0288e - 


- 10 


22.014084 


1.96165e 


- 10 


22.711070 


2.01973e 


- 10 


27.899896 


1.95265e 


- 10 


28.812177 


2.01571e 


- 10 


34.474785 


1.95008e 


- 10 


35.631088 


2.01904e 


- 10 


41.738699 


1.95577e 


- 10 


43.167733 


2.03125e 


- 10 


49.691673 


1.96711e 


- 10 


51.422281 


2.04763e 


- 10 


58.333258 


1.97959e 


- 10 


60.391382 


2.06229e 


- 10 


67.665387 


1.99038e 


- 10 


70.140636 


2.0725e - 


- 10 


77.957871 


1.99894e 


- 10 


79.828287 


2.0782e - 


- 10 


89.484519 


2.65662e 


- 10 


91.544035 


2.08206e 


- 10 



Table 3: The interlacing property of the first 15 purely imaginary eigenvalues for e = 0.48 and e = 0.5. 
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5 Conclusion 



We have proved that the operator L associated with the heat equation (jl.ip admits a closure in 
L 2 ([— 7r, 7r]) with a domain in H^ cv ([— n, it]) for |e| < 2. The spectrum of L consists of isolated eigenvalues 
of finite multiplicities. Furthermore, we have proved with the assistance of numerical computations 
that the set of eigenfunctions of the spectral problem (jl.3p for isolated eigenvalues is complete and 
minimal, but does not form a basis in H^ er ([— n, n]). By using the analytic function theory and the 
truncated matrix eigenvalue problem, we have approximated the eigenvalues numerically and showed 
that all eigenvalues of the spectral problem (|1.3|) are purely imaginary. 

We conclude therefore that the spectral problem (|1,3|) is not useful for direct analysis of well-posedness 
of the Cauchy problem for the heat equation (jl.ip . This is very unusual compared to the case of a 
standard heat equation 

h = —ho — efiQQ, t > 0, 



h(0) = h , (5 ' 1} 
subject to the periodic boundary conditions on 9 £ [—tt,tt]. By using the method of separation of 
variables and the completeness of the Fourier series in H^ er ([— it, 7r]), one can easily prove that there 
exists a unique solution of the heat equation (|5,ip for any ho G H^ er {[— it, it]) and any e G M. in the form 

h{t) = J2 Cne^e^-V, t > 0, (5.2) 

where the coefficients {c n } ne ^ are found in one and only one way from the initial data 

V/io G tfjUhvr, vr]) : h (9) = ]T c n e inx . (5.3) 

If e > 0, the Cauchy problem for the backward heat equation (|5.ip is ill-posed. If ho is only in 
Hp CT ({— 7r, 7r]), coefficients {c n } ne z decays to zero algebraically fast as n — > oo and the solution (15. 2p is 
singular for any t > 0. Therefore, the solution h(t) blows up in an infinitesimally small time. 

If e < 0, the Cauchy problem for the forward heat equation (|5.ip is well-posed. In this case, the 
solution h[t) becomes analytic in 9 G [— 7r, vr] for any t > even if ho is only in Hp eT ([—ir, 7r]). Therefore, 
there exists a constant C > such that ||^(i)||j?i er ([-7r,7r]) ^ C|l^o||i/i er ([-7r,7r]) for any t > 0. 

Unlike this classical situation, the series of eigenfunctions of the spectral problem (jl.3p are not 
applicable to construct a solution of the periodic heat equation (jl.ip . unless conditional convergence 
proposed in [4j can be adopted on a rigorous footing. Therefore, the PDE analysis of the Cauchy 
problem for the periodic heat equation (II. ip remains an open problem up to the date. 

Note: When the project was essentially complete, we became aware of a preprint [7], where similar 
results were obtained. In particular, the author of [7] proves that the spectral problem (jl.3p has no 
essential spectrum and, assisted with the numerical computations, illustrates that the eigenfunctions 
for all isolated eigenvalues do not form a basis. The analysis of [7| is based on the difference eigenvalue 
problem (|4.6p . which makes it different from our analysis. 
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A Spectrum of the linear operator L in weighted spaces 



The operator L can be rewritten in the Sturm-Liouville symmetric form 



cot|- 



1A L L - d 



tan — 
2 



sm( 



d9 



Let r{9) = | tan(0/2)| 1 / e be the weight of the Sturm-Liouville spectral problem 

-eL f(9) = Xr(9)f(9), 
acting on smooth functions f(9) on 9 £ [— 7r,7r] in the weighted space / £ L 2 ([—tt,tt]). 



(A.l) 



(A.2) 



Proposition 2 The operator Lq admits a self-adjoint extension in L 2 ([— 7r, it]) forO < e < 1, such that 
the spectrum of Lq is purely discrete, consists of a set of simple real eigenvalues {X n }nei with Ao = 0, 
X n = — X- n , Vn G N and lim A n = oo, and the eigenfunction f n {9) for X n > is identically zero on 

n— >oo 

6 G [-tt,0]. 



Proof. By Lemma[71 the eigenfunction f(6) of the spectral problem (|A.2|) for < e < 1 is in L^([— 7r,7r]) 
if and only if it is spanned by the fundamental solutions fi(0), /^(^) and f^iO), such that f(9) is 



bounded at 6 = and lim f(6) = 0. Let f(9) on ±9 £ [0,tt] be represented by 



icol ( - 



-i l/2e 



g±(x), cos 6 = x, 



where x £ [—1, 1]. Then, the spectral problem fjl .3j) is rewritten in the form 



dx 



(1 



4e(l - x 2 



Ag±(x) 

vr 



x^ 



(A.3) 



(A.4) 



If f{9) belongs to Lf([— tt, it]) for < e < 1, then g±(x) belong to L 2 p ([—1, 1]) with the weight function 
p[x) = -^t==5 . Since the symmetric spectral problem (|A.4p is self-adjoint in space L 2 p {[— 1,1]), its 
eigenvalues A are all real- valued. By Lemma [3] and the Sturm-Liouville theory [5], the eigenvalues 
{A n } ne z are all simple and symmetric with Ao = and X n = —X- n , Vn £ N, while the sequence 
{X n } n( zz is unbounded with lim A n = oo. By the standard Green's identity, 



x \\9+ 



2 



;i-x 2 )|<(x)| 2 dx + 



1 |ff + (x)| 2 dx 
x 4e(l -x 2 ) 



>0, 



(A.5) 



the eigenvalues X n of g+{x) are proved to be positive. If X n > is an eigenvalue for <?+(x), then 
— X n < can not be an eigenvalue for g~(x), such that g~{x) = on x £ [—1,1]. Therefore, f n {9) = 
on 9 £ [— 7T, 0] if A n > 0. By Lemma El the spectrum of L is purely discrete for |e| < 2, such that the 
spectrum of Lq in L 2 ([— 7r,7r]) is also purely discrete. □ 
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Remark 5 The operator Lq admits the same self-adjoint extension for any e > 1. However, this 
extension is not unique. Indeed, if e > 1, the eigenfunction f(9) of the spectral problem (|A,2p may 
exist in L^([— 7r, tt]) even if it is spanned by both fundamental solutions of Lemma [7] in the form 

f(9) = Ah{9) + Bf 2 (8) = A + f+(9) + B+f+(0) = A~^{9) + B~f-(9). 

It follows from the Sturm-Liouville problem ()A.2p that 



A 



r{9)\f{e)\ 2 d9 



tan — 
2 



sm{e)f{e)f{e)\\ 



+e £hG0 



sin(fl)|/'(0)| 2 ^. 



The second integral is finite for e > 1 only if i? = 0. In this case, the first term is computed explicitly 
for e > 1 as follows 

2 1 / 6 (A + B+ + A~B~) . (A.6) 

This term is zero if A + = A~ = for the self-adjoint extension described in Proposition [2l However, 
this choice is not unique, e.g. the alternative pairing B + = B~ = can also be applied. 



B Resonant poles of the Schrodinger operators 

Eigenvalues of the operator L coincide for |e| > | with resonant poles of the Schrodinger operators. 
To show this, we use the transformation (|3,5|) on the intervals ±9 £ [0, tt] and rewrite the spectral 
problem (jl.3p as the uncoupled spectral problems (|3.6p for f(9) = f±(t) on t S M. Let us apply the 
normalization condition /(tt) = /(— 7r) = 1 for the eigenfunction f(9) of the spectral problem (jl.3|) 
with A ^ R. Indeed, if either f(ir) = or /(— vr) = 0, then the identity (|A.5[) implies that A G R, 
but a real-valued A can not be an eigenvalue of the spectral problem (|1.3|) by Lemma [U Therefore, 
eigenfunctions f±(t) of the uncoupled problems (|3.6p satisfy the boundary conditions 

lim f±(t) = l, lim f ± (t)=a ± , (B.l) 

t—*— OO I— >oo 

where a 1 * 1 are uniquely defined. The function f{0) on £ [— tt,it} constructed from f±(t) on i S 1 is 
continuous at 9 = if a + = a~. 

Let e > ^ and define f±(t) = e^g±{t). The eigenfunctions g±(t) satisfy the linear Schrodinger 
equations 

{h ~ ed ') 9± = ±Asecht g ±> ( B - 2 ) 

but the boundary conditions (jB.ip show that g± ^ L 2 (M) and A is not an eigenvalue. In fact, the 
eigenfunctions g± (t) belong to the exponentially weighted L 2 -space, such that A is a resonance pole of 
the Schrodinger operators. Let us decompose the eigenfunctions by g±(t) = e~ ~^ +h±(t). By the theory 
of exponential asymptotics of solutions of the Schrodinger problems (|B.2|) . it follows that h± £ L 2 (R) 
if e > h. Let ho(t) = sechi G L 2 (R) and define the linear inhomogeneous problems for h±(t): 

Sf(X)h± = ±Xh (t), Sf(X) = ^ - e^ 2 ^ Asechi. (B.3) 

The operator S^(X) maps H 2 (R) to L 2 (M) and, if A ^ R, the kernel of Sf 1 is empty. The boundary 
condition a + = a~ is then equivalent to the zeros of the function 

H e (X) = X lim [(S+)- x (A) + (S^y^-X)] h {t). (B.4) 
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The function H £ (X) coincides (up to a multiplicative constant) with the analytic function F e (X) intro- 
duced in Corollary Q] for |e| < 2. Therefore, H t (X) represents a continuation of F e (X) from the domain 
|e| < 2 to the domain |e| > \. The function H e (X) is analytic in A G C and its roots give isolated 
eigenvalues of the spectral problem (|1.3|) with the account of their multiplicity. 

The function H e (X) can be simplified for A G £DL Let A = iuo G iM and represent h± = F(t) ± iG(t), 
where 

L e F = -ijj sechi G, L t G = oj secht F + uh (t), L e = — - ed 2 . (B.5) 
Therefore, we can define a real- valued function H e (uj) = H t (iuj) onwGf given by 

HJuj) = oj lim (L e + J 2 secht L' 1 secht)" 1 h (t). (B.6) 

t— >oo 

By multiplying the linear inhomogeneous equation 

(L e + oj 2 secht L^ 1 secht) G = ujho(t), 
by and integrating on i G M by parts, the function H t (uj) can be represented in the integral form 

/OO /'OO ^ 

sechtcft — uj 2 e^secht L~ l sechtG(t)dt. 
-00 J —00 

Let /io(*) = secht L7 1 sechte^. Then, 

/oo 
h (t)G(t)dt. 
-00 

This is yet another form of the real analytic function, zeros of which are equivalent to purely imaginary 
eigenvalues of the spectral problem (|1.3j) for |e| > \. 
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